Alternating-Spin Ladders in a Magnetic Field: New Magnetization Plateaux 
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We study numerically the formation of magnetization plateaux with the Lanczos method in 2-leg 
ladders with mixed spins of magnitudes (Si,S2) = (1,1/2) located at alternating positions along 
the ladder and with dimerization 7. For interchain coupling J' > and 7 = 0, we find normalized 
plateaux at m = 1/3 starting at zero field and m — 1 (saturation), while when 7 7^ is columnar, 
another extra plateau at m = 2/3 shows up. For J' < 0, when 7 < 7c(J') we find no plateau 
while for 7 > jc(J') we find four plateaux at m = 0, 1/3, 2/3, 1. We also apply several approximate 
analytical methods (Spin Wave Theory, Low-Energy Effective Hamiltonians and Bosonization) to 
understand these findings and to conjeture the behaviour of ferrimagnetic ladders with a bigger 
number of legs. 



PACS number: 76.50.+g, 75.50.Gg, 75.10.Jm 



I. INTRODUCTION 

Low dimensional spin systems has attracted many in- 
terests in recent years. It is the dominant effect of quan- 
tum fluctuation which distinguishes them from the clas- 
sical picture which is valid for higher dimensions. For 
instance it was argued by Haldanc jij that integer spin 
chains have an energy gap between the ground state and 
the first excited state. In the presence of magnetic field, 
the energy gap persists up to a critical field, equal to 
the gap. For 5=1 Anti- Ferromagnetic Heisenberg (AFH) 
chain the spectrum is gapless from the critical field up to 
the saturation field The 5=1/2 AFH chain remains 
gapless from zero field up to the saturation field, where 
the ground state is fully polarized |J. The situation 
might be more complex when one considers non-nearest 
neighbour interactions or non-homogeneous spin system 
by considering Bond- Alternations (BA) or different kind 
of spins. However a similar information may be obtained 
about the full spectrum by study the magnetization of 
system in the presence of a magnetic field. 

The study of the magnetization process in one- 
dimensional quantum spin systems has been spured since 
the discovery of quantized plateaux in the magnetization 
curves as a function of the applied mag netic field g, §, 
H . This behaviour is a sort of "magnetic analogy" of the 
Quantum Hall Effect for charge degrees of freedom and it 
has become an issue of both theoretical and experimen- 
tal interest. Likewise, ld-chains with alternating mixed 
spins have been shown to behave with unusual quantum 
properties: they show ferrimagnetic behaviour instead of 
antiferromagnetism |t]]- [jnj. For instance for low spin 
values of (S±,S2) ferrimagnetic chain (5i ^ S2) the low 
energy spectrum is gapless while there is a gap above this 
spectrum. Thus, a natural question arises as to whether 
ferrimagnetic systems are on equal footing with AF sys- 
tems as far as the formation of magnetization plateaux is 



concerned. Here a remark is in order. As an alternating 
spin chain has an ordered ground state (GS) with total 
spin 5 = -2"(5i — S2), then a natural plateau is expected 
even at zero field. Being this GS magnetization quite 
robust, a possible scenario for quantum ferrimagnets is 
that no other plateaux will appear, in addition to the 
saturation one. We shall show that this is not the case. 
Although there are recent studies of one-chain ferrimag- 
nets with magnetic fields E3] , jl3| , to our knowledge lit- 
tle is known about the onset of magnetization plateaux 
in ferrimagnetic ladders. A ferrimagnetic ladder may be 
regarded as the first step to consider the interaction be- 
tween ferrimagnetic chains as shown in Fig.l. This struc- 
ture exists in MnCu(pbaOH)(H20)-i where it is com- 
posed of parallel planes containing ferrimagnetic chains 
and in each plane the interaction between chains is as- 
sumed to be smaller than the intra-chain interactions 
fuj . In this paper we address the study of ferrimagnetic 
ladders under several situations of the Anti- (J' > 0) and 
Ferro-magnetic (J 7 < 0) interchain coupling and 7 dimer- 
ization patterns, which reveal a very reach structure of 
plateaux at fractional values of the full magnetization. 

We define the total magnetization M per unit length 
as follows, 

M=\j^Y,(S^*{n)+S^{n)), (1) 

a—l n—1 

where (n) denotes the quantum spin-1 at site n in 
the leg a = 1,2 of the ladder, similarly S^^n) for the 
spin-1/2 and the number of sites is N = 2 x L, L is 
the number of rungs. Alternatively we may use to, the 
normalized value of the magnetization with respect to 
the saturation value as to = M/M sat . We want to study 
this observable in a 2-leg ferrimagnetic ladder with the 
following Heisenberg Hamiltonian (i = 2n + a), 
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+ [1 + 7 (Q) (* + 1)] S^ (i + l)-S^ } (z + 2) ) 

1,2 L 

+ Jl E E s i Q) (^) • s 2%) - ^ otai , (2) 

n=l 

where ft is the magnetic field, J and J' are coupling 
constants along the legs and the rungs respectively, and 
(n) is the dimerization patterns. We consider two dif- 
ferent dimerization patterns : i) CBA (Columnar Bond 
Alternation, Fig|l]-(a)) for which ^ a \n) = (-1)" +1 7 
PI, and ii) SBA (Staggered Bond Alternation, Fig.0- 
(b)) where ^ a \n) = (-l)«+ 1 +« 7 ||. We use periodic 
boundary conditions along the legs of the ladder. 

The paper is organized as follows. In the next section 
we show the application of spin wave theory to derive the 
saturation field. In Sec. Ill we obtain a low-energy ef- 
fective Hamiltonian based on Quantum Renormalization 
Group (QRG) formalism to describe the magnetization 
curve and give a conjecture for a 2n-leg ferrimagnetic lad- 
der. We present an Abelian Bosonization in Sec. IV to 
support our conjecture and express the necessary condi- 
tion for the appearance of the magnetization plateaux in 
a 2rt-leg ferrimagnetic ladder. In Sec. V we will confirm 
the results of the analytical methods using a numerical 
method (Lanczos method) applied to a 2-leg ladder un- 
der a variety of situations. Finally, we summerize our 
results as a conclusion in Sec. VI. 



II. SPIN WAVE THEORY 

We start considering the case J' > 0. Then, since 
the GS is an ordered state with non-vanishing magne- 
tization, it is natural to apply SWT to study the low 
energy properties of these ladders and predict correctly 
their ferrimagnetic behaviour. Now we extend this anal- 
ysis including the effect of the external magnetic field. 
For weak magnetic fields we expect that the ferrimag- 
netic GS will be little perturbed so that it is still a good 
choice to make a SWT analysis around a partially po- 
larized state with magnetization Sq S = ^{S\ — S%). At 
h = the 2-leg ladder presents two pairs of branches, say 
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,(±) 



(fc; h = O),^^ h = 0), in the dispersion relation, 



namely: 



\k; 0) = i(d - C 2 + yl{C 1 +C 2 ) 2 -lCl {± \k)), 
<4 ±} (*; 0) = \{C2 - C x + y/(C 1 + C 2 ) 2 -4C 2 3 i± \k)), (3) 

d = (2 J + J')Si ; C 2 = (2 J + J')S 2 , (4) 
C { 3 ±} (k) = ^S!S 2 (2J cos k± J') . (5) 
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FIG. 1. A picture of the dimerization patterns considered 
in a 2-leg mixed spin ladder: a) Columnar Bond Alternation 
(CBA), b) Staggered Bond Alternation (SBA). 

For J' = the two chains decouple and the branches 
collapse by pairs reproducing the one-chain result. When 
h is applied the spectrum (Bl) is modified as follows: 
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(6) 



After doing this analysis we find that the Bogoliubov 
transformation to diagonalize the SWT Hamiltonian is 
independent of h. This fact implies that we cannot see 
the plateaux of the curves M = M(h) using this SWT. 

For strong magnetic fields the previous ferrimagnetic 
GS is not a good choice to perturb around. It is more 
natural to choose a fully polarized (ferromagnetic) state 
with total magnetization Sq S = ^{S\ + S 2 ). We find 
that the SWT analysis yields the following spectrum: 

tt^; 0) = \{~{Ci + C 2 ) ± ^(C 1+ C2) 2 -^+)(fc)), 

4 ±} (*; 0) = + Ca) ± V / (Ci+C 2 ) 2 -£>(")(*:)), 

D (±) (k) = 16JS 1 S 2 (J± Jcosfc + J')(l Tcosfc). (7) 

The inclusion of the magnetic field merely shifts this spec- 
trum by h: 



Sl^ (k; h) = Qf (Jfc; h = 0) + h , 
n { s\k; h) = n ( g \k; h = 0) + h 



(±)/ 



(8) 



Again, we are not able to detect magnetization plateaux. 
However, now we can give a reliable prediction for the 
location of the saturation field h sat . This happens when 



(-) 



(k = 0; h sat ) = 0, giving us the following relation for 
1, S 2 = 1/2 and normalizing J = 1: 



\{2 + J'). 



(9) 
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We shall find that this condition is verified numerically. 
For decoupled chains J' = we reproduce the one-chain 
results [OJ. Similarly, for ferromagnetic coupling J' < 
we find the following saturation field, 
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FIG. 2. Magnetization plateaux for alternating 2-leg lad- 
der: 7 = 0, J' = 0.2, 1.0, 2.0. There are two plateaux at 
m = 1/3 , 1. 



III. LOW ENERGY EFFECTIVE HAMILTONIAN 

In order to explain qualitatively the onset of magneti- 
zation plateaux in a 2-leg ferrimagnetic ladder we can use 
an RG transformation to map the original model (||) onto 
a one-chain model with known properties. To perform 
the block RG transformation we choose as the building 
block the unit cell made up of 2 spins (Si = 1, S2 = 1/2). 
The block Hamiltonian is then h B = J'S v S 2 -h{S(+S%). 
For simplicity, we assume a strong coupling regime J' 3> 
J. This allows us to place h B on the vertical rungs of 
the 2-leg ladder so that the whole block Hamiltonian is 
the rung Hamiltonian (H B — H r ) and the inter-block 
Hamiltonian is the leg Hamiltonian (H BB = Hi). We 
need to distinguish 2 regimes of magnetic field to con- 
struct the truncation operator, a) For low fields, namely 
h < 3/2J', the lowest lying states of hP is a spin-1/2 
doublet. Then we use this to map (||) onto the following 
low-energy effective Hamiltonian (J = 1): 

L L 

H eff = ~LJ' - § 2_j 7n S n ' S n+1 ~ h S n > 
ri—l n— 1 

2 

7; = E( 1+ ^ (Q) M)> ( n ) 

a=l 
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FIG. 3. Magnetization curves for SBA configuration with 
J' = 2.0, 7 = 1.0,0.5 which shows only two plateaux at 
m = 1/3 , 1. 

H e ff is a spin (S' =) 1/2 Ferromagnetic Heisenberg chain 
in the presence of an external magnetic field, which shows 
a ferromagnetic phase with M = 1/2 (m = 1/3). 
b) For high fields, namely h > 3/2 J, the lowest two 
energy states of h B have now spin quantum numbers 
I St = 1/2, = +1/2), |3/2,+3/2). Upon projection 
of (H) onto this subspace we arrive at the following effec- 
tive Hamiltonian (J = 1): 

Heff = -(^ + h)L+^J2 < + Hxxz+h , (12) 

n— 1 

where Hxxz+h is an AF spin 1/2 Heisenberg model in 
a magnetic field with coupling parameters given by: 

2 L 

HxXZ+h = - ^2 T4( ^2 S 'n S n+l + ^ S 'n S 'n+l) ~ h 'n S n 
n—1 fi—x,y 




In the absence of dimerization 7 = or when the dimer- 
ization pattern is SBA, we have 7^ = 2 constant and 
thus the effective field hi is also uniform. In this case we 
can use the knowledge of the Hxxz+h model to predict 
the existence of 2 plateaux. It is known that this model 
has a critical value for hi where a phase transition occurs 
between the Partially Magnetized (PM) to the Saturated 
Ferromagnetic (SF) states. This critical line in the A — h' 
phase diagram is given by |l6| h' c = A + 1 (for hi > 0). 
Firstly, we find that the relation between the magneti- 
zation M of the original ferrimagnetic ladder and the 
magnetization Mxxz of the effective XXZ model is: 

M = l + M xxz , (14) 

Now it is possible to derive the dependence of M on h. By 
using the relation in (113) we find two critical fields given 
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by h c 



l c l — f + f J' an( i ^c2 = if i . 

then M = 1/2, for h cl < h < h c2 then 1/2 < M < 3/2 
and for h > h c2 then M = 3/2. Therefore, we find a good 
qualitative picture of the M (/i)-plot with two plateaux 
at normalized values of m = 1/3 and m = 1 |l7j . This 
agrees with the numerical results in Fig. 1 and Fig. 2 
using Lanczos. 

Although the knowledge of known models obtained in 
Eqs.( pd| , ^3| ) was used to explain the magnetization curve 
one may however continue using QRG for Eqs.([ll]|l3|) to 
obtain the RG-flow of coupling constants and finally the 
magnetization curve. In the case of Eq.(p"l]) the Hamilto- 
nian will be projected to another ferromagnetic Heisen- 
berg model with larger spins. By repeating this proce- 
dure one will arrive at the classical ferromagnetic Heisen- 
berg fixed point. In the other case in Eq. ([l3|) , the RG- 
flow will describe both the partially magnetized and the 
fully polarized state of the XXZ+h model [jl8). Needless 
to say that the final result will not change. 

Furthermore, if we restrict ourselves to a strong cou- 
pling analysis, we can extend this QRG analysis in a 
simple way to ferrimagnetic ladders with a number of 
legs given by n; = 2n in order to preserve the structure 
of (Si, S f 2)-cells in both directions. Then, in each rung 
the highest total spin is n(Si + S 2 ) and when we ap- 
ply a h field we find that the allowed values of M are 
-n(Si + S 2 ), —n(Si + S 2 ) + 1, • • • , n(Si + S 2 ) - 1, n(Si + 
S 2 ). Upon normalization to the highest value, we get 
that the allowed values of m are: 



I J' such that for h < h cl and Right fermionic fields V>l iTO 2 -) 



1,-1 
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1 



n(Si+S 2 )'' 



n(Si + S 2 )' 



+1. (15) 



Now, if we particularize to the case Si = 1, S 2 — 1/2, we 
find -1, -1 + 2/3n, . . . , +1 - 2/3rc, +1. We can readily 
check that in this set the values 1/3 and 1 are always 
present. Moreover, we may conjeture the existence of 
additional plateaux if the number of legs increases. 



= 1,2 such that 
^,W,ci ^ e «^W' Ci (2;) + e- ikFX ^f ,Ci {x). Then, the 
bosonization proceeds in a standard fashion representing 
these fermions by vertex operators in terms of chiral bo- 



son fields ^l'r as V'l 'r{ x ) — e Tl ^ L - R respectively, 
and R is the compactification radius of the boson. In this 
low energy analysis, all potentially massive modes will 
decouple and only the massless bosons will contribute to 
this effective description. This is the case of the total 

diagonal boson $ = £, i=12 Y^Ui^L + <^fl)- In 
terms of this boson, the translational symmetry by one 
unit cell of the original lattice chain is realized as a total 
shift $ -> $ + 2(Si + S 2 )k F 2irR. On the other hand, the 
Fermi momentum is modified by a shift of the magne- 
tization as hp = §(1 — m) (the magnetic field acts as a 
chemical potential.) With these two conditions, a leading 
operator of the form cos($/i?) will be compatible only if 
(Si + S 2 )(l — m) = integer. This is the application of the 
OYA condition to ferrimagnetic systems. 
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IV. ABELIAN BOSONIZATION 

We may re-obtain this quantization condition using a 
heuristic argument based on Abelian Bosonization, as in 
the original OYA argument This has been carried 
out with full detail for standard AF spin- 1/2 ladders in 
refs. fl9fl . For simplicity, let us consider first a one-chain 
ferrimagnet (J' > 0) with spins (Si, £2) i n the unit cell 
and no dimerization 7 = 0. This unit cell will play the 
role of an effective site in the bosonization. Firstly, the 
spin chain is fermionized by introducing colored fermions 

(i) c- 

xpn'' % i = 1, 2, with color index Cj = 1, . . . , 2Si to de- 
scribe the spin degrees of freedom of Si and S2, respec- 
tively. We make the usual assumption about the con- 
tinuum limit: only excitations close to the Fermi points 
are important. Thus, we introduce slowly varying Left 



FIG. 4. Magnetization plateaux for CBA configuration 
with J' = 2.0, 7 = 0.8, 0.5. A new plateau appears at 
m = 2/3, its width increases by increasing 7. 

In our case, we have Si = 1,S 2 = 1/2 and this yields 
only m = 1/3, 1. This agrees with the QRG analysis Jl~7| ] 
and with numerics ]l3[ |. We may attempt to go beyond 
this one-chain case to a set of n 2-leg alternating ladders 
coupled together. In each 2-leg ladder we assume a strong 
coupling limit (J' 3> 1) and place the unit cell (Si,S2) 
on the rungs as in the RG analysis. Then, we assume 
a weak coupling (J' 1) among these set of n ladders. 
Repeating the argument yields the following quantization 
condition for magnetization plateaux: 



T u (Si + S 2 ) n (1 - m) = integer 



(16) 
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where we have introduced a "translational unit" T u to ac- 
count for the translational symmetry of the ground state. 
In the original Hamiltonian we have T u — 2 but due to 
non-perturbative effects it can be different in the real 
ground state This condition ([l6| ) (with T u = 1) is 
precisely equivalent to eq. (|T^) JL9|. We remark that 



this condition (16) appears as a necessary condition, but 
it is not sufficient. Additionally, the compactification ra- 
dius should be greater than its critical value. Thus, we 
consider this quantization condition as a conjeture and 
we have applied the Lanczos method to test it. 



NUMERICAL RESULTS: 
METHOD 



THE LANCZOS 



Now we present our numerical studies of 2-leg ladders 
with a number of sites N — 2 x L where L = 4, 6, 8, 10 
(due to the constraint of periodic BC's). We apply Lanc- 
zos to each 5 z -sector of the Hilbert space from S z = 
to S z — y(£>i -I- S%). We always set J — 1 and make 
J' range both in AF and Ferromagnetic regimes. The 
dimerization parameter 7 is also varied for both CBA 
and SBA patterns. 

Firstly, we consider AF couplings (J 1 > 0) and vanish- 
ing dimerization (7 = 0). In Fig.H we show 3 plots of 
the magnetization curve for weak J' = 0.2, intermedi- 
ate J' = 1.0 and strong coupling J' = 2.0. For weak 
coupling we find two plateaux at m = 1/3, 1 and values 
of h c \, h C 2 = h sat very close to the pure Id ferrimagnet. 
This amounts to a smooth evolution from the one-chain 
case. Moreover, at strong and intermediate couplings, 
the same picture holds true which supports our effective 
Hamiltonians ([TT|) , (p~2|) based on the low-energy analysis. 
We have also checked numerically that the value of the 
saturation field h sat fits nicely the SWT result in eq. (|^), 
as clearly shown in FigJ|. Another remarkable aspect of 
our plots is that they show a square-root behaviour of the 
magnetization m— m c ~ \Jh? — h% in the transient region 
near the critical points, as is also the case for standard 
AF spin- 1/2 ladders Next, as a check of our numer- 
ical method, we consider SBA dimerization with 7=1 
and J' — 2.0. In this situation the 2-leg ladder degener- 
ates onto a one-chain ferrimagnet with coupling J = 2.0 
and in Fig.|| we reproduce exactly the existence of two 
plateaux and the value of the saturation field h sat = 3 
(in units J = 2.) We find that this SBA picture is robust 
by plotting the case 7 = 0.5. 

In order to see the effect of different dimerization patterns 
in the formation of the magnetization plateaux, we have 
also studied the CBA dimerization configuration. We 
plot the case CBA for J' = 2.0 and 7 = 0.5, 0.8 in Fig]|. 
Interestingly enough, we clearly find a new plateau at 
to = 2/3 in addition to the previous ones at to = 1/3 and 
to = 1. This agrees with eq. ( |l6| ) with T u = 2. The width 
of this new plateau grows with increasing dimerization. 



The new plateau shows the appearance of an energy gap 
between the two energy bands in the spectrum of CBA 
configuration. This picture is true for CBA configuration 
(J' > 0) when 7^1. By setting 7=1 the ladder will be 
decomposed to isolated 4-sites plaquette which does not 
reflect the properties of a system in the thermodynamic 
limit (N — > 00). 

Secondly, we consider Ferromagnetic couplings (J' < 0) 
with different dimerizations. In this case the CBA pat- 
tern is more interesting for we have found that it exhibits 
two phases, one gapless for 7 < j c and another gapped 
for 7 > 7c, in the J' — 7 diagram with a phase-separating 
line which goes as j c ~ 0.42 when J' — » — 00 fl5 |. 

In Fig.|| we consider a plot for J' = —1.0 and 7 = 0.2 
which happens to be in the gapless phase. Accordingly 
with this behaviour, there is no plateau at to = and the 
transient region immediately starts at h = 0. Moreover, 
the previously robust 1/3-magnetization plateau is now 
absent. The reason for this is because now the GS has no 
longer total spin ^(Si —52) but 0. This is completely in 
agreement with the effective Hamiltonian which has been 
obtained for this region in jni) . The effective Hamiltonian 
for CBA configuration (J' < 0) and 7 < 7 C is a bond- 
alternation spin 3/2 AFH model in the gapless phase [pl| . 
It has been shown that the gapless phase of spin 3/2 
AFH chain is in the universality class of a spin 1/2 AFH 
chain ]2C| ]. As it has been noted in the introduction the 
model has no plateau in the magnetization curve up to 
the saturation field. 

Finally, we consider the opposite case of parameters 
corresponding to the gapped phase, namely J' = —1.0 
and 7 = 0.8, in Fig.|. 
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FIG. 5. Magnetization curves for ferro CBA ladders with 
J' = —1.0, 7 = 0.2. There is no plateau which reveals no gap 
in the full spectrum. 
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This plot shows a very rich structure. It shows four 
magnetization plateaux at m = 0,1/3,2/3,1. Since the 
model has an energy gap in this region the m = plateau 
extends until the magnetic field closes the gap in the 
spectrum h c — A. More interesting is now the m — 
1/3 plateau, which is not natural for now the GS is not 
ferrimagnetic when J' < 0. In fact, this plateau does 
not start at h = as in Figs|^-[|. To explain the nature 
of these plateaux we can resort to a strong coupling RG 
analysis and map the original Hamiltonian (|^) onto a 
AFH one-chain of spins 3/2 with dimerization (see eq.(5) 
in jl!|) and apply the OYA argument. We have also 
checked that the saturation field h sat is also correctly 
predicted by SWT in eq. ©. 

The situation of SBA at J' < is similar to Fig.|. This 
can be understood by the effective low-energy Hamilto- 
nian which is a homogenous spin 3/2 AFH chain [fL5| . 
The model is gapless for the whole range of magnetic 
field which results to no plateaux. 




FIG. 6. Magnetization plateaux for ferro CBA configura- 
tions with J' — —1.0, 7 = 0.8. Four plateaux appear in this 
region at m = , 1/3 , 2/3 , 1. 



VI. CONCLUSIONS 

We have addressed the study of alternating (Si,^)- 
spin ladders in the presence of an external uniform mag- 
netic field in a fixed direction. The low energy prop- 
erties of this problem is relevant in understanding the 
behaviour of ferrimagnetic systems in magnetic fields as 
compared with the properties exhibited by more standard 
antiferromagnetic ladders. 

In summary, we have shown that mixed-spin systems 
show extra magnetization plateaux in addition to the nat- 



ural ones at m = 1/3 and 1 and we have presented a 
complete account for the diversity of plateaux-patterns 
when we vary the interchain and dimerization couplings. 

We have undertaken a complete numerical study us- 
ing the Lanczos method of the 2-leg ladders with alter- 
nating spins Si = 1/2, S2 = 1, for several regimes of 
the interchain coupling constant J', both antifcrro- and 
ferromagnetic, and for several types of dimerization pat- 
terns (CBA and SBA). In the antiferromagnetic regime 
(J' > 0) the SBA configuration shows only the natural 
plateaux at m = 1/3 and 1 while the CBA configura- 
tion shows an extra plateau at m = 2/3. In the ferro- 
magnetic regime (J' < 0) the SBA configuration shows 
no plateau where the model remains gapless from the 
zero field up to the saturation field. The CBA configura- 
tion in the J' < regime has more interesting features. 
For 7 < 7 C the model is gapless for all range of mag- 
netic field and results to no plateau, while for 7 > j c the 
model is gapful. Then the first plateau appears at m = 0. 
There are also two other plateaux at m — 1/3 and 2/3 
which are expressed by the dimerized S—3/2 AFH ef- 
fective model. The trivial plateau at m = 1 counts the 
number of plateaux to 4 in this case. 

The magnetization curves obtained in this way provide 
a whole picture of the response of mixed-spin systems to 
external magnetic fields. 

We have complemented these numerical studies with 
additional analytical methods such as Spin Wave The- 
ory (in two fashions), effective low-energy mappings of 
the type used in Quantum RG methods and Abelian 
Bosonization. Although none of these analytical methods 
provide a complete understanding of the problem, they 
however give us some useful partial knowledge that fits 
nicely into our findings obtained with numerical meth- 
ods. 

Concerning the possible experimental realization of the 
magnetic ferriladders studied in this work, we would like 
to point out that the field of organic chemistry is a likely 
candidate for synthesizing this type of materials [^l] ^3| . 
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